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I give a brief discussion of possible sources of high energy neutrinos of astrophysical origin over the energy range 
from ~ 10 12 eV to ~ 10 25 eV. In particular I shall review predictions of the diffuse neutrino intensity. Neutrinos 
from interactions of galactic cosmic rays with interstellar matter are guaranteed, and the intensity can be reliably 
predicted to within a factor of 2. Somewhat less certain are intensities in the same energy range from cosmic rays 
escaping from normal galaxies or active galactic nuclei (AGN) and interacting with intracluster gas. At higher 
energies, neutrinos will definitely be produced by interactions of extragalactic cosmic rays with the microwave 
background. With the discovery that gamma ray bursts (GRB) are extragalactic, and therefore probably the most 
energetic phenomena in the Universe, it seems likely that they will be copious sources of high energy neutrinos. 
Other sources, such as AGN and topological defects, are more speculative. However, searches for neutrinos from 
all of these potential sources should be made because their detection would have important implications for high 
energy astrophysics and cosmology. 



1. INTRODUCTION 

The technique for constructing a large area (in 
excess of 10 4 m 2 ) neutrino telescope has been 
known for more than two decades |lj|. The pio- 
neering work of the DUMAND Collaboration led 
to the development of techniques to instrument a 
large volume of water in a deep ocean trench with 
strings of photomultipliers to detect Cherenkov 
light from neutrino- induced muons (2) . Locations 
deep in the ocean shield the detectors from cosmic 
ray muons. The second generation of high energy 
neutrino telescope such as AMANDA || located 
deep in the polar ice cap at the South Pole, and 
NT 200 in operation in Lake Baikal, Siberia [Q, 
have demonstrated the feasibility of construct- 
ing large area experiments for high energy neu- 
trino astronomy. The next generation telescopes, 
such as the planned extension of AMANDA, ICE- 
CUBE @, and ANTARES §, may have effective 
areas of 0.1 km 3 , or larger, and be sufficiently 
sensitive to detect bursts of neutrinos from ex- 
tragalactic objects and to map out the spectrum 
of the diffuse high energy neutrino background. 
In this paper I focus on possible astrophysical 
sources of neutrinos contributing to the diffuse 
high energy neutrino background from ~ 1 TeV 
to the GUT scale. 



2. COSMIC RAY INTERACTIONS 
WITH MATTER 

There will definitely exist a diffuse galactic 
neutrino background due to interactions of the 
galactic cosmic rays with interstellar matter. 
The spectrum of cosmic rays is reasonably well 
known, as is the matter distribution in our galaxy. 
Estimates of the neutrino intensity have been 
made by Silberberg and Shapiro 0, Stecker ||, 
Domokos et al. ||, Berezinsky et al. fjfjfl , and In- 
gelman and Thunman |Tl)] , and the more recent 
predictions are shown in Fig. [lj. The differences 
of about a factor of 2 between the predictions are 
accountable in terms of the slightly different mod- 
els of the interstellar matter density, and cosmic 
ray spectrum and composition used. Also shown 
is the atmospheric neutrino background as esti- 
mated by Lipari In addition, there will be 
a very uncertain background (not plotted) due to 
charm production (see refs. p3| , [l4| for a survey 
of predictions). 

Somewhat less certain is the flux of neutrinos 
from clusters of galaxies. This is produced by pp 
interactions of high energy cosmic rays with intr- 
acluster gas. Berezinsky et al. [|l5| have made 
predictions of this, and I show in Fig. ^| their 
estimates of the diffuse neutrino intensity due 
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Figure 1. Neutrinos from cosmic ray interactions with the interstellar medium (upper curves for t = 

0°, 6 = 0°, lower curves for b = 90°): - - Domokos et al. |J; Berezinsky et al. 

- Ingelman and Thunman pJj| . The band with vertical hatching shows the range of atmosheric 
neutrino background |l2| as the zenith angle changes from 90° (highest) to 0° (lowest). Neutrinos from 
cosmic ray interactions with the microwave background: — • — Protheroe and Johnson |^3|| for 

£ max = 3 x 10 20 eV and 3 x 10 21 eV; Hill and Schramm ||]; assuming the 

highest energy cosmic rays are due to GRB according to Lee |25j . 



to interactions of cosmic rays produced by nor- 
mal galaxies and AGN together with an upper 
limit based on assuming the observed 7-ray back- 
ground results from 7r° production. Later esti- 
mates by Colafrancesco and Blasi 16 are also 
shown. 



3. COSMIC RAY INTERACTIONS 
WITH RADIATION 

Moving to higher energies, cosmic rays above 
~ 10 20 eV will interact with photons of the cos- 
mic microwave background radiation (CMBR) 
|l7| , |l8|| . Again, we know that both ingredients 
exist (the highest energy cosmic ray detected has 
an energy of 3 x 10 20 eV fl9|] , and at least 6 cos- 
mic rays have been detected above 10 20 eV by the 
AG AS A array [^0) ) , and so pion photoproduction 
at these energies will occur, resulting in a diffuse 
neutrino background (Stecker B). However, the 



intensity in this case is model-dependent because 
it is not certain precisely what the origin of the 
highest energy cosmic rays is, and whether in fact 
they are extragalactic, although this seems very 
probable (see for a discussion of the high- 
est energy cosmic rays). One of the most likely 
explanations of the highest energy cosmic rays 
is acceleration in Fanaroff- Riley Class II radio 
galaxies as suggested by Rachen and Biermann 
Protheroe and Johnson Q have repeated 
Rachen and Bicrmann's calculation in order to 
calculate the flux of diffuse neutrinos and 7-rays 
which would accompany the UHE cosmic rays, 
and their result has been added to Fig. |l[ Any 
model in which the cosmic rays above 10 20 eV are 
of extragalactic origin will predict a high energy 
diffuse neutrino intensity probably within an or- 
der of magnitude of this at 10 19 eV. For example, 
I show an earlier estimate by Hill and Schramm 
p4[ . Also shown is an estimate by Lee [E5| of the 
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Figure 2. Neutrinos from cosmic ray interac- 
tions in clusters of galaxies (Berezinsky et al. 
|15[): lower hatched area - normal galaxies; up- 
per hatched area - AGN; dashed line - upper 
bound from 7-ray data. Large hatched region - 
comic ray interactions with intergalactic medium 
in clusters of galaxies jl6| 



diffuse neutrino intensity estimated in a model in 
which the highest energy cosmic rays have their 
origin in sources of gamma ray bursts. 

4. GAMMA RAY BURSTS 

Gamma ray bursts (GRB) are observed to have 
non-thermal spectra with photon energies extend- 
ing to MeV energies and above. Recent identi- 
fication of GRB with galaxies at large redshifts 
(e.g. GRB 971214 at z = 3.42 |f|) show that 
the energy output in 7-rays alone from these ob- 
jects can be as high as 3 x 10 53 erg if the emis- 
sion is isotropic, making these the most energetic 
events in the Universe. GRB 980425 has been 
identified with an unusual supernova in ESO 184- 
G82 at a redshift of z = 0.0085 implying an en- 
ergy output of 10 52 erg (^tJ. These high energy 
outputs, combined with the short duration and 
rapid variability on time-scales of milliseconds, re- 



quire highly relativistic motion to allow the MeV 
photons to escape without severe photon-photon 
pair production losses. The energy sources of 
GRB may be neutron star mergers with neutron 
stars or with black holes, collapsars associated 
with supernova explosions of very massive stars, 
hyper-accreting black holes, hypernovae, etc. (see 
|||,|2j| for references to these models). 

The relativistic fireball model of GRB 
(Meszaros and Rees |Q) provides the framework 
for estimation of neutrino fluxes from GRB. A 
relativistic fireball sweeps up mass and magnetic 
field, and electrons are energized by shock acceler- 
ation and produce the MeV 7-rays by synchrotron 
radiation. Protons will also be accelerated, and 
may interact with the MeV 7-rays producing neu- 
trinos via pion photoproduction and subsequent 
decay at energies above ~ 10 14 eV pl| , p3[ . Accel- 
eration of protons may also take place to energies 
above 10 19 eV, producing a burst of neutrinos at 
these energies by the same process |34[] . These en- 
ergetic protons may escape from the host galaxy 
to become the highest energy cosmic rays J|5],[36| . 
Additional neutrinos due to interactions of the 
highest energy cosmic rays with the CMBR will 
be produced as discussed in the previous section. 

For a sufficiently intense GRB, it may be pos- 
sible to identify neutrinos from individual GRB. 
Integrating over all GRB in the Universe, Wax- 
man and Bahcall |u],[52| have predicted the dif- 
fuse neutrino intensity, and this has been plotted 
in Fig. ^ with a steepening at 10 16 eV, and with 
a continuation to higher energies as suggested by 
Vietri §§. 

5. ACTIVE GALACTIC NUCLEI 

The 2nd EGRET catalog of high-energy 7-ray 
sources jSTJ contains over 40 high confidence iden- 
tifications of AGN, and all appear to be blazars 
(radio-loud AGN having emission from a rela- 
tivistic jet closely aligned to our line of sight). 
Since the publication of the 2nd catalog, the num- 
ber of blazars detected by EGRET has increased 
to nearly 70 (see refs. f|^|3{| for reviews). TeV 
emission has been observed from three blazars, 
the BL Lac Objects Mrk 421, Mrk 501 and 1ES 
2344+514 jf0|. Clearly, the 7-ray emission is as- 
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Figure 3. 
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Diffuse neutrinos from 
solid line Waxman & Bahcall pl| ] 
sion to higher energies suggested by Vietri p4| 
(thick chain line) . Diffuse pr/ (lower) and pp + pry 
(upper) neutrinos from blazars: - — Mannheim 
Q (jry) Model A (lower curve), {pj+pp) Model 
B (upper curve); - - Protheroe Q x25% 

- see text; pry Halzen & Zas [pil. 



sociated with AGN jets. Blazars appear also to 
be able to explain about 25% of the diffuse 7-ray 
emission |4l| |, and models where 7-ray emission 
does not originate in the jet are unlikely to con- 
tribute significantly to the diffuse 7-ray (and neu- 
trino) intensity (see Protheroe and Szabo |l2| and 
references therein for predictions for non-blazar 
AGN). Several of the EGRET AGN show 7-ray 
variability with time scales of ~ 1 day |Q at GeV 
energies, and variability on time scales of ~ 1 
hour or less 1 44 ij| has been observed at ~ 1 TeV 
for some BL Lacs. These variability timescales 
place important constraints on the models, and 
not all models developed so far are consistent 
with this. I shall survey the neutrino emission 



predicted in blazar models irrespective of this, as- 
suming they may be made to accommodate the 
latest variability measurements. 



Most theoretical work on 7-ray emission in 
AGN jets involved electron acceleration and in- 
verse Compton scattering, and these models will 
predict no neutrinos. In proton blazar models, 
protons are accelerated instead of, or as well as, 
electrons. In this case interactions of protons with 
matter or radiation would lead to neutrino pro- 
duction. In some of the proton blazar models 
energetic protons interact with radiation via pion 
photoproduction (see e.g. for references and 
a discussion of pry interactions). This radiation 
may be reprocessed or direct accretion disk radi- 
ation p6| , or may be produced locally, for exam- 
ple, by synchrotron radiation by electrons accel- 
erated along with the protons [E7 48 1. Pair syn- 



chrotron cascades initiated by photons and elec- 
trons resulting from pion decay give rise to the 
emerging spectra, and this also leads to quite ac- 
ceptable fits to the observed spectra. These mod- 
els can produce neutrinos and also higher energy 
radiation than electron models because protons 
have a much lower synchrotron energy loss rate 
than electrons for a given magnetic environment. 
In both classes of model, shock acceleration has 
been suggested as the likely acceleration mecha- 
nism (see ^jj for references). 

By appropriately integrating over redshift and 
luminosity in an expanding universe, using a lu- 
minosity function (number density of objects per 
unit of luminosity) appropriate to blazars, and us- 
ing the proton blazar models to model the 7 ray 
and neutrino spectra one can estimate the dif- 
fuse neutrino background expected from blazars. 
In Fig. H I have added intensities of (i/f, + v^) 
predicted in proton blazar models by Mannheim 
@, Protheroe @ (x0.25 as only ~25% of 7-ray 
background is due to AGN [^4j - original calcula- 
tion assumed 100%) and Halzen and Zas |l9|]. For 
some of these models expected muon rates have 
been calculated WtWi. 



6. TOPOLOGICAL DEFECTS 

Finally, I discuss perhaps the most uncertain 
of the components of the diffuse high energy neu- 
trino background, that due to topological defects 
(TD). In a series of papers ]53|45^], TD have 
been suggested as an alternative explanation of 
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Figure 4. The result of || for M x c 2 = 10 
GeV, a magnetic field of 10 -9 gauss, and p = 2, 
normalized the spectrum of "observable particles" 
(nucleons, photons, electrons) to the 3x 10 11 GeV 
data H]). H = 75 km s" 1 Mpc" 1 and q = 0.5 
are assumed. 



the highest energy cosmic rays. In this scenario, 
the observed cosmic rays are a result of top-down 
cascading, from somewhat below (depending on 
theory) the GUT scale energy of - 10 16 GeV 
|57[ , down to 10 11 GeV and lower energies. These 
models put out much of the energy in a very flat 
spectrum of neutrinos, photons and electrons ex- 
tending up to the mass of the "X-particles" emit- 
ted. 

Protheroe and Stanev |5^] argue that these 
models appear to be ruled out by the GeV 7- 
ray intensity produced in cascades initiated by 
X-particle decay for GUT scale X-particle masses. 
The 7-rays result primarily from synchrotron ra- 
diation of cascade electrons in the extragalactic 
magnetic field. Fig. ||, taken from rcf. |5q| , shows 
the neutrino emission for a set of TD model pa- 
rameters just ruled out according to Protheroe 
and Stanev |5^] for a magnetic field of 10~ 9 G 
and X-particle mass of 1.3 x 10 14 GeV. Clearly for 
such magnetic fields and higher X-particle masses 
(e.g. GUT scale), TD cannot explain the highest 
energy cosmic rays. Indeed there is evidence to 
suggest that magnetic fields between galaxies in 



clusters could be as high as 10~ 6 G [|55|. However, 
for lower magnetic fields and/or lower X-particle 
masses the TD models might explain the highest 
energy cosmic rays without exceeding the GeV 
7-ray limit. For example, Sigl et al. |60| show 



that a TD origin is not ruled out if the extra- 
galactic field is as low as 10~ 12 G, and Birkel & 
Sarkar adopt an X-particle mass of 10 12 GeV. 
Yoshida et al. |32| investigate various TD scenar- 
ios with GUT scale masses, and their predicted 
neutrino fluxes are generally higher than those of 
Sigl et al. but such high neutrino intensi- 

ties are likely to be excluded because the 7-rays, 
due to cascading even in a 10~ 12 G field, would 
probably exceed the GeV flux. The intensities are 
compared in Fi g. |e| . A novel feature of the work of 
Yoshida et al. 62 is the inclusion of interactions 
of high energy neutrinos with the 1.9 K cosmic 
neutrino background, and this can be important 
at the very highest energies. 

I emphasize that the predictions summarized 
in Fig. H are not absolute predictions, but the 
intensity of 7-rays and nucleons in the resulting 
cascade is normalized in some way to the high- 
est energy cosmic ray data. It is my opinion that 
GUT scale TD models are neither necessary nor 
able to explain the highest energy cosmic rays 
without violating the GeV 7-ray flux observed. 
The predicted neutrino intensities are therefore 
extremely uncertain. Nevertheless, it is impor- 
tant to search for such emission because, if it is 
found, it would overturn our current thinking on 
the origin of the highest energy cosmic rays and, 
perhaps more importantly, our understanding of 
the Universe itself. 



7. DISCUSSION 

Very recently, Waxman and Bahcall |52| have 
used some arguments based on the observed cos- 
mic ray spectrum to obtain an upper bound to 
high energy neutrinos from astrophysical sources. 
Their argument hinges on sources of astrophysi- 
cal neutrinos being sources of the highest energy 
cosmic rays which happen also to produce neutri- 
nos by P7 interactions. Hence, except for sources 
with a very high optical depth for protons, the 
maximum neutrino intensity will be about 10% of 
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Figure 5. Neutrinos from topological defects: 
Stanev || (10 141 GeV, 10"" 
(2 x 10*GeV, 10" 12 gauss); 
(10 16 GeV, gauss). 
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— TD model just ruled out according to Protheroe & 
gauss); - - TD model just allowed according to Sigl et al. [6C] 
Birkel & Sarkar Q (10 12 GeV, gauss); Yoshida et al.[62] 



the extragalactic (~ E~ 2 ) component of the high- 
est energy cosmic rays. Examining AGN models, 
they find that predictions for proton blazar mod- 
els exceed their bound. 

In the case of AGN, they also suggest that the 
optical depth to pj at ~ 10 19 eV must be much 
less than 1 (to enable TeV 7-rays to escape with- 
out significant 77 pair production losses), with 
the consequence that the ultra high energy cosmic 
ray production far exceeds the ultra high energy 
neutrino production, pushing the neutrino upper 
bound even lower. 

TeV 7-rays, however, have so far only been seen 
from 3 blazars, and it is by no means certain that 
TeV 7-rays are emitted by all blazars and so high 
P7 optical depths are not necessarily ruled out 
(note, however, that the infrared background lim- 
its how far away one can observe objects at TeV 
energies t33|). Also, in at least one of the proton 
blazar models [H the optical depth of protons to 
P7 at 10 19 eV is high because the proton direc- 
tions are isotropic in the jet frame whereas the 
radiation field is highly anisotropic, coming from 



near the base of the jet, and the photons cascade 
down to TeV energies where the 77 optical depth 
along the jet direction is low because of the radi- 
ation being anisotropic. Admitedly, neutrons are 
produced in a fraction of P7 interactions, and the 
neutrons escape as cosmic rays, and so the effec- 
tive optical depth for nucleons can not exceed ~ 1 
by much, and so it is probable that this proton 
blazar model is ruled out. 

The main argument relating the neutrino up- 
per bound to the observed ultra high energy cos- 
mic ray flux relies on the cosmic rays of energy 
10 19 eV being able to reach Earth from AGN dur- 
ing the Hubble time. There is evidence to suggest 
that magnetic fields between galaxies in cores of 
clusters (the most likely place to find an an AGN) 
could be as high as 10~ 6 G |9). With such high 
magnetic fields it is not obvious that 10 19 eV pro- 
tons will reach us from most AGN contributing 
high energy neutrinos. 

Thus, I believe that the "upper bound" is 
model dependent, and that its calculation is com- 
plicated by cosmic ray propagation effects. While 
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Figure 6. Grand Unified Neutrino Spectrum - a 
personal opinion about the predicted neutrino in- 
tensities: thick solid lines - certain; long dashed 
lines - almost certain; short dashed lines - spec- 
ulative; dotted lines - highly speculative. 



I would certainly classify the higher AGN fluxes 
as speculative, or highly speculative, I believe the 
lower ones are not ruled out by the argument of 
Waxman and Bahcall. Nevertheless, the work of 
Waxman and Bahcall is very important in re- 
minding us that for any model used to predict 
high energy neutrino fluxes we must check that it 
does not overproduce cosmic rays. 

Plotting a representative sample of the diffuse 
flux predictions from Figs. [|, [2| || and |s| in the 
same figure one has a "grand unified neutrino 
spectrum" (with apologies to Ressell and Turner 
64 1) . This is shown in Fig. ^| where I have labelled 
the various curves as "speculative" , "highly spec- 
ulative", "certain" or "almost certain". These 
labels reflect my own personal opinion or preju- 
dice and should not be taken too seriously - other 
opinions are equally valid. 



8. PROSPECTS FOR OBSERVATION 

With the construction in the relatively short 
term of 0.1 km 2 neutrino telescopes, and in the 
longer term of 1 km 2 detectors, it is useful to esti- 
mate the signals expected due to various possible 
neutrino intensities. At high energies, electron 



neutrinos may also be detected through the re- 
sulting cascade, and this is particularly important 
when looking for horizontal air showers, for ex- 
ample with the proposed AUGER detectors fl62f . 
Several estimates of event rates have been made 
for various energy thresholds, or for horizontal air 
showers due to neutrino interactions (including v e 
and D e , see e.g. Gandhi et al. |52[). 

To illustrate how the (y^ + v^) signals ex- 
pected from different astrophysical neutrino spec- 
tra would be detected by telescopes with different 
energy thresholds, I have made approximate esti- 
mates of the event rates as a function of minimum 
muon energy using the P 1/ ^ fl (E u , E™ ln ) function 
given in Fig. 2 of ref. @ for E™ in = 1 GeV, mod- 
ified for other £™ m values in a way consistent 
with that given for E^ in = 1 TeV. The effects 
of shadowing for vertically upward-going neutri- 
nos have been included using the shadow factor 
S(E V ) given in Fig. 20 of ref. |51j| . I have esti- 
mated the expected neutrino induced muon signal 
for four representative neutrino intensities. The 
vertically upward-going and horizontal muon sig- 
nals are shown separately for each case in Fig. |^ 
together with the atmospheric neutrino induced 
muon signals for the two directions. As can be 
seen, the highest signals would be due to the pro- 
ton blazar models, with several events per year 
expected in a 0.1 km 2 detector. However, one 
should be cautious as these intensities are some- 
what speculative (as discussed earlier) . Detection 
of muon signals in one year from the other inten- 
sities estimated would be marginal for a 0.1 km 2 
detector, but achievable with a 1 km 2 detector. 
Detection of transient neutrino signals, correlated 
with observations of the same source in photons 
(e.g. GRB, AGN) should therefore be the goal 
of high energy neutrino astronomy in the short 
term. 

One should consider the consequences for as- 
trophysical neutrinos of the discovery of the os- 
cillation of atmospheric v u , probably into u T , 
by Super-Kamiokande J65[ with an oscillation 
length of A osc ~ 10 3 (£'/GeV) km. On an as- 
trophysical scale, the oscillation length A osc ~ 
3 x 10 -11 (EyTeV) kpc is very small, and integrat- 
ing contributions to the neutrino intensity over 
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Figure 7. Atmospheric neutrino induced muon 
signal (thick solid lines) and putative astrophys- 
ical neutrino induced muon signals expected for 
the following (y^ + v^) intensities: GRB inten- 
sity |3l| extended to 10 20 eV (dotted curves); 



P7 proton blazar |48| (short dashed curves); TD 
model just ruled out according to ref. [[58| (solid 
curves); interactions of the highest energy cosmic 
rays with the CMBR §§ for E max = 3 x 10 21 eV 
(long dashed curves). Upper curves show horizon- 
tal signals, lower curves show vertical (upward) 
signals. 



astrophysical dimensions one would naively ex- 
pect the (Vfi + Dfj.) flux to be 50% lower (assuming 
sin 2 26 = 1), and to be accompanied by a similar 
{v T + i> T ) flux. The unique signature for detection 
of tau neutrinos has been discussed in ref. [ ]66[ . 

Neutrino astronomy is developing during an era 
in which exciting discoveries are being made in 
other areas of high energy astrophysics. These 
include detection of rapidly varying TeV 7-ray 
signals from AGN, discovery that GRB are extra- 
galactic and probably the most energetic phenom- 
ena ocurring in the Universe today, and detection 
at Earth of cosmic rays with energies well above 
10 20 eV opening the question of whether their 
origin is through particle acceleration at radio 
galaxies or GRB, or from topological defects left 
over from the big bang. Hadronic processes may 
have a role in all these phenomena, and search- 



ing for high energy neutrinos may lead to greater 
understanding of the highest energy phenomena 
in the Universe. Clearly, for this to happen in 
parallel with the other observations rapid devel- 
opment of neutrino telescopes with sensitive areas 
of ~ 1 km 2 or larger, operating over a wide range 
neutrino energies, is essential. 
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